Search for: All records

Abstract Very long baseline interferometry observations reveal that relativistic jets like the one in M87 have a limb-brightened, double-edged structure. Analytic and numerical models struggle to reproduce this limb-brightening. We propose a model in which we invoke anisotropy in the distribution function of synchrotron-emitting nonthermal electrons such that electron velocities are preferentially directed parallel to magnetic field lines, as suggested by recent particle-in-cell simulations of electron acceleration and the effects of synchrotron cooling. We assume that the energy injected into nonthermal electrons is proportional to the jet Poynting flux, and we account for synchrotron cooling via a broken power-law energy distribution. We implement our emission model in both general relativistic magnetohydrodynamic (GRMHD) simulations and axisymmetric force-free electrodynamic (GRFFE) jet models and produce simulated jet images at multiple scales and frequencies using polarized general relativistic radiative transfer. We find that the synchrotron emission is concentrated parallel to the local helical magnetic field and that this feature produces limb-brightened jet images on scales ranging from tens of microarcseconds to hundreds of milliarcseconds in M87. We present theoretical predictions for horizon-scale M87 jet images at 230 and 345 GHz that can be tested with next-generation instruments. Due to the scale-invariance of the GRMHD and GRFFE models, our emission prescription can be applied to other targets and serve as a foundation for a unified description of limb-brightened synchrotron images of extragalactic jets. 

Abstract We study the spherical accretion of magnetized plasma with low angular momentum onto a supermassive black hole, utilizing global general relativistic magnetohydrodynamic simulations. Black hole-driven feedback in the form of magnetic eruptions and jets triggers magnetized turbulence in the surrounding medium. We find that when the Bondi radius exceeds a certain value relative to the black hole’s gravitational radius, this turbulence restricts the subsequent inflow of magnetic flux, strongly suppressing the strength of the jet. Consequently, magnetically arrested disks and powerful jets are not a generic outcome of the accretion of magnetized plasma, even if there is an abundance of magnetic flux available in the system. However, if there is significant angular momentum in the inflowing gas, the eruption-driven turbulence is suppressed (sheared out), allowing for the presence of a powerful jet. Both the initially rotating and nonrotating flows go through periods of low and high gas angular momentum, showing that the angular momentum content of the inflowing gas is not just a feature of the ambient medium, but is strongly modified by the eruption and jet-driven black hole feedback. In the lower-angular-momentum states, our results predict that there should be dynamically strong magnetic fields on horizon scales, but no powerful jet; this state may be consistent with Sgr A* in the Galactic center. 

Abstract Massive elliptical galaxies harbor large amounts of hot gas (T≳ 10⁶K) in their interstellar medium (ISM) but are typically quiescent in star formation. The jets of active galactic nuclei (AGNs) and Type Ia supernovae (SNe Ia) inject energy into the ISM, which offsets its radiative losses and keeps it hot. SNe Ia deposit their energy locally within the galaxy compared to the larger few ×10 kiloparsec-scale AGN jets. In this study, we perform high-resolution (512³) hydrodynamic simulations of a local (1 kpc³) density-stratified patch of the ISM of massive galaxies. We include radiative cooling and shell-averaged volume heating, as well as randomly exploding SN Ia. We study the effect of different fractions of supernova (SN) heating (with respect to the net cooling rate), different initial ISM density/entropy (which controls the growth timet_tiof the thermal instability), and different degrees of stratification (which affect the freefall timet_ff). We find that SNe Ia drive predominantly compressive turbulence in the ISM with a velocity dispersion ofσ_vup to 40 km s⁻¹and logarithmic density dispersion ofσ_s∼ 0.2–0.4. These fluctuations trigger multiphase condensation in regions of the ISM, where $\min \left(t_{\mathrm{ti}}\right)/t_{\mathrm{ff}}\lesssim 0.6\exp \left(6{\sigma }_s\right)$, in agreement with theoretical expectations that large density fluctuations efficiently trigger multiphase gas formation. Since the SN Ia rate is not self-adjusting, when the net cooling drops below the net heating rate, SNe Ia drive a hot wind which sweeps out most of the mass in our local model. Global simulations are required to assess the ultimate fate of this gas. 

Abstract Low-collisionality plasma in a magnetic field generically develops anisotropy in its distribution function with respect to the magnetic field direction. Motivated by the application to radiation from accretion flows and jets, we explore the effect of temperature anisotropy on synchrotron emission. We derive analytically and provide numerical fits for the polarized synchrotron emission and absorption coefficients for a relativistic bi-Maxwellian plasma (we do not consider Faraday conversion/rotation). Temperature anisotropy can significantly change how the synchrotron emission and absorption coefficients depend on observing angle with respect to the magnetic field. The emitted linear polarization fraction does not depend strongly on anisotropy, while the emitted circular polarization does. We apply our results to black hole imaging of Sgr A* and M87* by ray tracing a GRMHD simulation and assuming that the plasma temperature anisotropy is set by the thresholds of kinetic-scale anisotropy-driven instabilities. We find that the azimuthal asymmetry of the 230 GHz images can change by up to a factor of 3, accentuating (T_⊥>T_∥) or counteracting (T_⊥<T_∥) the image asymmetry produced by Doppler beaming. This can change the physical inferences from observations relative to models with an isotropic distribution function, e.g., by allowing for larger inclination between the line of sight and spin direction in Sgr A*. The observed image diameter and the size of the black hole shadow can also vary significantly due to plasma temperature anisotropy. We describe how the anisotropy of the plasma can affect future multifrequency and photon ring observations. We also calculate kinetic anisotropy-driven instabilities (mirror, whistler, and firehose) for relativistically hot plasmas. 

ABSTRACT Recent observations with JWST have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the big bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox$$^{\it HR}$$, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments (FIRE) project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox$$^{\it HR}$$ re-simulates the cosmic volume ($L=22.1$ cMpc) of the original FIREbox run with eight times higher mass resolution ($$m_{\rm b}\sim {}7800\, M_\odot$$), but with identical physics, down to $$z\sim {}6$$. FIREbox$$^{\it HR}$$ predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at $$z\sim {}6{\!-\!}14$$, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox$$^{\it HR}$$, the SFE–halo mass relation for intermediate mass haloes ($$M_{\rm halo}\sim {}10^9{\!-\!}10^{11}\, {\rm M}_\odot$$) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE–halo mass relation lead to a larger contribution from lower mass haloes at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE–halo mass relation inferred from FIREbox$$^{\it HR}$$ allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at $$z\gt 12$$ will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn. 

ABSTRACT In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. Active galactic nucleus jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work concludes from simulations of a halo with $$10^{14} \,\mathrm{ M}_\odot$$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $$10^{12}-10^{15}\, {\rm M}_{\odot }$$ haloes using high-resolution, non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without violating observational constraints. On the contrary, if we scale the energy flux based on the total cooling rate within the cooling radius, strong interstellar medium cooling dominates this scaling, resulting in a jet flux exceeding what is needed. Among the three jet types, the CR-dominant jet is most effective in suppressing cooling flows across all surveyed halo masses due to enhanced CR pressure support. We confirm that the criteria for a successful jet model work across a wider range, encompassing halo masses of $$10^{12}-10^{15} {\rm M_\odot }$$. 

Abstract In 1977, Blandford and Znajek showed that the electromagnetic field surrounding a rotating black hole can harvest its spin energy and use it to power a collimated astrophysical jet, such as the one launched from the center of the elliptical galaxy M87. Today, interferometric observations with the Event Horizon Telescope (EHT) are delivering high-resolution, event-horizon-scale, polarimetric images of the supermassive black hole M87* at the jet launching point. These polarimetric images offer an unprecedented window into the electromagnetic field structure around a black hole. In this paper, we show that a simple polarimetric observable—the phase ∠β₂of the second azimuthal Fourier mode of the linear polarization in a near-horizon image—depends on the sign of the electromagnetic energy flux and therefore provides a direct probe of black hole energy extraction. In Boyer–Lindquist coordinates, the Poynting flux for axisymmetric electromagnetic fields is proportional to the productB^ϕB^r. The phase ∠β₂likewise depends on the ratioB^ϕ/B^r, thereby enabling an observer to determine the direction of electromagnetic energy flow in the near-horizon environment experimentally. Data from the 2017 EHT observations of M87* are consistent with electromagnetic energy outflow. Currently envisioned multifrequency observations of M87* will achieve higher dynamic range and angular resolution, and hence deliver measurements of ∠β₂closer to the event horizon as well as better constraints on Faraday rotation. Such observations will enable a definitive test for energy extraction from the black hole M87*.